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A global, six-level primitive equation model was numeri- 
cally integrated to examine the interaction of a planetary 
wave with a growing mid-latitude baroclinic wave. The 
purpose of this study was to determine the contribution of 
the baroclinic wave to the initiation of a monsoon surge. 
Three experiments were run. Experiment I depicted the 
effect of the planetary wave on the southward migration of 
cold air. The cyclogenesis of a mid-latitude baroclinic 
wave was shown in Experiment II. Experiment III combined 
the growing baroclinic wave with the stable planetary wave 
in an attempt to simulate a cold surge. Characteristics of 
a monsoon surge were present in the numerical simulation but 
not to the extent or detail of observed surges. The results 
indicated the important role of planetary wave advections 
in the formation of cold surges. 
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I. 



INTRODUCTION 



The Asian winter monsoon has been an area of very active 
research over the past few years. Winter monsoons are 
dramatic manifestations of land-sea-air interaction. These 
phenomena consist of a mixture of all scales of motions, and 
the various interactions are still not well understood. 
Although the fundamental structure of the monsoon has 
frequently been observed, the dynamical mechanisms are 
still unknown. Recent areas of investigation include the 
role of the winter monsoon in mid-latitude tropical inter- 
action, monsoonal variation and the forcing mechanisms for 
small scale monsoonal variations. 

The monsoon is a three-dimensional, planetary scale wind 
regime which exhibits a strong seasonal dependence. It 
exhibits two distinct seasonal low-level wind flow patterns- 
a winter outflow from a cold continental anticyclone and a 
summer inflow into a continental heat low. The winter anti- 
cyclone and the summer heat low result from land-sea tempera 
ture differences and the varying position of the sun. The 
shapes of the continents and their topographies also play 
a role in producing the considerable regional and temporal 
variability of the monsoons. 

Though centered far to the north, the Siberian anti- 
cyclone dominates the Southeast Asian winter (Ramage, 1971). 
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The low heat capacity of the land relative to the open ocean 
to the east and south ensures the development of the cold 
anticyclone with surface flow from land to ocean. A com- 
pensating flow from ocean to land occurs in the upper 
troposphere. The Himalayas block the southward movement of 
the extremely cold air from the Siberian anticyclone. As the 
anticyclone becomes established over Siberia, a surface 
cyclone develops over Indonesia and Malaysia as a major 
planetary scale convective area is shifted eastward from its 
summer position near India (Baker, 1983). 

Chang and Lau (1980) analyzed the monthly mean velocity 
potential field at 200 mb for December 1974 and showed a 
broad region of divergence over Indonesia. This broad area 
of divergence is indicative of deep rising motions and large 
scale overturning in the atmosphere. A region of maximum 
convergence over northern China, in the vicinity of the 
Siberian anticyclone, was also noted. 

A strong Hadley cell dominates the winter tropospheric 
circulation of this region. This thermally direct circu- 
lation involves the descent of cold air over the Siberian 
anticyclone, the southward flow of air at the surface, the 
ascent of warm air over Indonesia and a northward return 
flow at 200 mb. 

Although the large scale features of the monsoon are 
identifiable from year to year, there is large interannual 
variability in strength and location. Velocity, temperature 
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and pressure fields vary on short time scales. One such 
variation is the monsoon surge, also known as the "cold 
surge." These events occur at intervals of several days 
to two weeks during the winter monsoon (Lim and Chang, 1981). 

Observational studies by Danielson and Ho (1969), Chang 
et al . (1979), Murakami (1979), Chang and Lau (1980,1982) 

and Chang et al. (1983) discuss the characteristics and 
effects of the monsoon surge at the synoptic and planetary 
scales. The characteristics and effects described in these 
studies are summarized in the following paragraphs. 

Over northern Asia there is intense, deep cold air 
advection associated with the passage of a mid-tropospheric 
trough. The trough intensifies rapidly over northern Japan. 
The baroclinicity increases with the anticyclogenesis over 
Siberia and cyclogenesis in the strong baroclinic zone off 
the east China coast. As the anticyclone intensifies, 
subsidence increases which strengthens the east Asian Hadley 
circulation. This intensification increases the upper-level 
ageostrophic flow, causing an acceleration of the east Asian 
subtropical jet near Japan. 

The initiation of the monsoon surge occurs when the upper- 
level trough moves eastward and the pressure gradient along 
the east China coast increases. When the pressure gradient 
force exceeds the Coriolis force, the normal geostrophic 
balance is broken, resulting in cross-isobaric flow from 
high to low pressure. This accelerating northeasterly 
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low-level flow toward lower pressure is the monsoon surge. 

The surge reaches the equatorial South China Sea in 12-24 
h. In the equatorial regions, convection intensifies within 
one to two days, sustaining the already enhanced Hadley 
circulation. The upper-level outflow from the convection 
accelerates eastward and westward along the equator as well 
as northward. The surge ends when the mid-latitude trough/ 
ridge pattern moves far enough east, thus diminishing the 
strong subsidence over the anticyclone. 

The monsoon surge has two stages separated by a few 
hours to one day depending on the location of the observing 
site (Lim and Chang, 1981; Chang et a]^. , 1983) . The pres- 
surge surge is the first stage. It is the leading edge of 
the downward accelerating air with a propagation speed 
characteristic of internal gravity waves. A frontal passage 
is the second stage which moves with advective speeds and 
is defined by a sharp decrease in surface dewpoint. 

The cause or nature of monsoon surge forcing has been 
examined by Lim and Chang (1981) and Baker (1983). Linearized 
shallow-water equations on an equatorial beta-plane were 
used by Lim and Chang (1981) in an analytical study of 
monsoon response to simulated anticyclogenesis. They sug- 
gested that the large features of synoptic scale motions 
in the east Asian winter monsoon could be explained by 
pressure forcing using equatorial beta-plane dynamics. No 
consideration was given to orography, a variable mean zonal 
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wind, the vertical structure of monsoon surges or the upper 
tropospheric response to the surge. 

Baker (1983) used a global, six-level primitive equation 
model to examine the interaction of a mid-latitude baro- 
clinic wave with topography. A well-developed baroclinic 
wave could initiate a monsoon surge, however, the surges 
were weak with limited extent. Baker (1983) concluded that 
other forcing mechanisms were required to simulate the 
surges observed during the winter monsoon. 

Shaffer et a]^. (19 84) have shown that cold surge events 

are frequently preceded by rapid growth in the planetary 
waves. It appears that the amplified planetary wave helps 
advect the cold air toward the equator. The purpose of this 
study is to examine the interaction of a planetary wave 
with a growing baroclinic wave. The simulations will be 
carried out with a global, six-level primitive equation 
model. The same baroclinic jet will be used, as was em- 
ployed by Baker (1983). The baroclinic wave, which is unstable 
with the jet, will have wave number 8. The planetary wave 
will be represented by wave number 4, which is stable with 
the jet. 
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II . 



DESCRIPTION OF MODEL 



A version of the UCLA potential enstrophy conserving 
general circulation model is used in these experiments. The 
atmosphere is adiabatic, dry and frictionless. The governing 
equations of the model are the primitive equations (listed 
in Appendix A) . The model design and dynamic processes 
are described by Arakawa and Lamb (1977,1981) . The version 
used for these experiments is summarized below. 

A. THE HORIZONTAL GRID 

Arakawa ' s scheme C is used. The prognostic variables-- 
wind (u,v) , potential temperature (o) and terrain pressure 
(“)--are staggered horizontally as shown in Figure 2.1. This 




Figure 2.1. Horizontal distribution of variables in 
the UCLA model grid (Arakawa and Lamb, 
1977) . 



grid arrangement requires more computer time than an unstag- 
gered grid, but it provides a much better structure for the 
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shorter waves arising from geostrophic adjustment, smaller 
scale heating, nonlinear interaction or topographic effects. 
Grid resolution is 2.5° longitude by 4.0° latitude. A 90° 
sector of the globe extending from pole to pole is used as 
the computational domain. Cyclic continuity is imposed at 
the east and west boundaries. 

Convergence of the meridians toward the poles normally 
requires a very short time step to ensure computational 
stability. Instead of shortening the time step, longitu- 
dinal averaging is done on the terms governing gravity wave 
propagation (Arakawa and Lamb, 1977). The averaging is done 
by expanding the zonal pressure gradient and zonal mass flux 
into a Fourier series around each latitude circle. The 
amplitude of each wave component is then reduced by the 
factor S. The gradients of the meteorological fields are 
smoothed instead of their Fourier expansions. Application 
of the smoothing operator, S, occurs at high latitudes where 
S < 1. The smoothing operator is defined as: 

Q = AA cos ({) 

A(|) sin(mAA/2) 

(A list of symbol definitions is contained in Appendix B.) 

B. THE VERTICAL GRID 

The model equations are written using the vertical sigma 
coordinate : 
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a 



p -Pt 

— n — 



( 2 . 2 ) 



where : 



n = Ps - Pf 

is the surface pressure, and p^ is the pressure at the 
top of the model atmosphere (p^ = 200 mb) . The variables 
are vertically staggered in six equally spaced sigma layers 
(Figure 2.2). The surface and terrain pressure are defined 



p = 200mb — — a = 0 

u.v.e 



u.v.G 



u V 6 



U.V.0 



— U.V0 



p = Pj a = 1 



6-=0 

a 

& 

& 

& 

cr 

cr=0 



Figure 2.2. Vertical distribution of variables in a 

six-layer sigma coordinate system. Solid 
lines are o-levels; dashed lines represent 
the middle of the layer (Arakawa and 
Lamb, 1977) . 



on the lower boundary, o = 1. The vertical velocity, o, 
vanishes at the upper and lower boundaries and is diagnostically 
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derived on sigma levels. Wind (u,v) and potential tempera- 
ture (0) are computed at mid- layer. 

C. FINITE DIFFERENCING SCHEMES 

The model uses two time differencing methods--an Euler- 
backward scheme and the centered time differencing (or 
leapfrog) scheme. The model integration proceeds in a 
series of one Eurler-backward time step followed by four 
centered time steps. The Euler-backward scheme eliminates 
the computational mode of the centered scheme and selectively 
damps high frequency waves. The time step for these model 
integrations is 360 s. 
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III. MODEL INITIALIZATION 



An analytical initialization scheme developed by Hayes 
(1985) is used. In this scheme, the initial mean wind is 
specified analytically using the Bickley jet described by 
Haltiner and Williams (1980) . The initial surface pressure 
is derived from the wind using the geostrophic wind equation, 
and the initial temperature is derived using the thermal 
wind equation. In this research, a stable planetary wave 
and a synoptic scale perturbation are superimposed on a 
horizontally and vertically varying mean current. 

A. INITIALIZATION OF THE MEAN FLOW 

A baroclinic jet is located at 45° N and has a maximum 
velocity of approximately 46 m/s at 200 mb. The initial 
mean surface flow is given by the following equations: 



It has the same latitudinal variation as the jet with a 
maximum speed of 5 m/s at 45° N decreasing to zero at the 
poles. The v-component of the wind is initially set to 
zero. The mean upper-level wind is given by: 



Us(<|)) 




(3.1) 



Vs (4)) 



0 . 



(3.2) 
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(3.3) 



f P ( o , (}) ) ) 



sech^ (Y ) ( 



^n(Pg/p) 






v^(<l>/P(<^/4)) ) 



0 



(3.4) 



The wind varies in the vertical as the natural logarithm of 
pressure . 

The mean surface pressure is calculated by integrating 
the geostrophic wind equation: 



where p^ ((J)^ = 45° N) = 1013.25 mb and Tg(<}>) represents the 

mean surface temperature. Because the latitude structure of 

T (0) is unknown, p (())) is computed iteratively until the 
s s 

adjustment of T to the NACA standard for that p is less 

s s 

than 0.01° K. 

The mean surface temperature is calculated by integrating 
the geostrophic thermal wind equation which gives: 



The integrals are computed using the Simpson Rule approxima- 
tion. T(p,(})q) is obtained from the NACA standard atmosphere: 



Pg (0) = exp [Hn Pg (4>o) 



(a/R) / fUg ((J))/Tg (4))d())] 



(3.5) 



T(p,())) 
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t(Pq,<})q) - rz 



(3.7) 



where : 




FR/g 



Z 



The value of Tq , 288° K, is taken fromtheNACA standard 
atmosphere using a standard surface pressure of 1013.25 mb. 
(Constants are defined in Appendix B.) 

B. INITIALIZATION OF THE HAURWITZ WAVE 

The initialization of the Haurwitz wave is accomplished 
in a subroutine developed by Dr. T. Rosmond.^ The analyti- 
cal solution for a Rossby wave on the surface of a sphere 
is described by Haurwitz (1940). 

The Haurwitz wave has a zonal wave nurobef 4 (Figure 
3.1) . The maximum meridional velocity is 16 m/s. The u and 
V components are given by the following equations : 



V' = (3 

a cos 4) 8 X 

The equations of motion are used to derive the perturbation 
pressure : 



^Personal communication, present affiliation Naval 
Environmental Prediction Research Facility, Monterey, CA. 
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Figure 3.1. Sea- level pressure field for a stationary 
Haurwitz wave (3.10). Maximum meridional 
velocity is 16 m/s. Contour interval is 
20 mb. 
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p' 



C2 (u)+J^) cos ( Bt+mA) sin (f>p^ 



(3.10) 



n-m+1 m . 
n(n+l) ^n+l"^ 



where p^(sin (|)) is the associated Legendre polynomial and w 
is the mean flow. 

The equations listed above are for a Haurwitz wave 
superimposed on a mean flow. The Bickley jet and the 
Haurwitz wave schemes both have mean flow contributions, and 
one of the mean flows needs to be removed. Instead of using 
the formal analytical approach, the mean u, v and p for 
each J (latitude strip) is subtracted from the Haurwitz solu- 
tion to provide the Haurwitz wave contribution. This method 
allows the sea-level pressures to be computed at the same 
time as the winds for the Haurwitz wave and the Bickley jet. 
Otherwise, the winds would have to be computed and then the 
pressures. The error in this approach should not have a 
significant effect on the results. 

The Haurwitz wave u, v and p contributions are added to 
the Bickley jet and mean surface flow values. The surface 
values are spread to the upper levels using the equations 
described in the previous section. 

C. THE PERTURBATION 

The perturbation is a barotropic disturbance consisting 
of a longitudinal sine wave which is used for baroclinic 
instability studies. The maximum latitudinal amplitude in 
the meridional direction occurs at 45° N. The following 
equation describes the perturbation: 
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(3.11) 



(A sin (mA) sin ^(2cf))) 

The hydrostatic equation and the geostrophic wind 
equation are used to derive; 



' 

RT 



pf 



RT 



^ (A sin (mA) sin^ ( 2(J)) ) 



(3.12) 



u 



1 3 $’ 



af ^ 3(1) 



- —(A sin (mA) 2 sin (4(1)) ) 

O. 



(3.13) 



V' 



1 3 $ 



afp^ cos (}) 3(1> a cos (p 



^ (m Acos (mA) sin^ (2(1>) ) (3.14) 



where A = a*l m/s. These perturbation quantities are added 
to the mean current and integrated forward in time to 
produce the forecast. 
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IV. ANALYSIS OF RESULTS 



The hypothesis is that a cold surge results from the 
interaction of a growing wavenumber 8 cyclone moving east- 
ward and a quasi-stationary stable forced wave number 4. 

Hayes (1985) showed that the basic baroclinic current, which 
is given in Equations (3.1) and (3.3), is unstable for wave 
number 8 but is stable for wave number 4. The stationary 
wave could be partially forced by the Himalayas . Three 
experiments were run to verify the hypothesis. Experiment 
I shows the effect of the stable wavenumber 4 on the south- 
ward migration of cold air. The effect of cyclogenesis in 
wave, number 8 is depicted Experiment II. Experiment III 
is a combination of Experiments I and II. 

A. EXPERIMENT I 

In this experiment, the model is initialized with the 
stable wave number 4 and the baroclinic mean flow described 
in Chapter III. 

The baroclinic jet is confined to the northern hemisphere 
to minimize forcing from the southern hemisphere. The 
model is integrated forward in time to 60 h. Sea-level 
pressure fields for 12, 24, 36, 48 and 60 h are shown in 
Figures 4.1 to 4.5. The potential temperature fields at 
a = 11/12 for the same times are displayed in Figures 4.6 to 
4.10. The corresponding wind vectors are shown in Figures 
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Figure 4.1. Sea-level pressure field at 12 h for 

Experiment I. Contour interval is 5 mb. 
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Figure 4.2. Sea-level pressure field at 24 h for 

Experiment I. Contour interval is 5 mb. 
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Figure 4.3. Sea-level pressure field at 36 h for 

Experiment I. Contour interval is 5 mb. 
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0 




Figure 4.4. Sea-level pressure field at 48 h for 

Experiment I. Contour interval is 5 mb. 
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Figure 4.5. Sea-level pressure field at 60 h for 

Experiment I. Contour interval is 5 mb. 
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Figure 4.6. Potential temperature field at 12 h for 
Experiment I. Contour interval is 5° K. 
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Figure 4.7. Potential temperature field at 24 h for 
Experiment I. Contour interval is 5° K. 
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Figure 4.8. Potential temperature field at 36 h for 
Experiment I. Contour interval is 5° K. 
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Figure 4.9. Potential temperature field at 48 h for 
Experiment I. Contour interval is 5° K. 
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Figure 4.10. Potential temperature field at 60 h for 
Experiment I. Contour interval is 5° K. 
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4.11 to 4.15. The 90° sectors are plotted twice to produce 
hemispheric displays. 

In the first twelve hours, the effects of the baroclinic 
jet on the wave are readily apparent (Figure 4.1). The wave 
has moved westward approximately 5° compared to 10° in the 
southern hemisphere. The westward movement is caused by the 
beta effect. A slight NE-SW tilt in the sea-level pressure 
field is apparent. The absence of a corresponding tilt 
in the southern hemisphere is consistent with the absence of 
the jet. In the northern hemisphere, the tilt is in the 
same sense as the shear on the south side of the jet indi- 
cating that the wave will lose energy to the mean flow. The 
wave should damp. 

The pressure gradient has weakened as the cyclone and 
anticyclone are reduced in amplitude by 5 mb. However, at 
the same time, the amplitude of wave number 4 in the poten- 
tial temperature field (Figure 4.6) has increased. The 
305° K isotherm has moved south to near 10° N. The jet has 
advected that field eastward to a position 45° out of phase 
with the pressure field. Northeasterly winds of approxi- 
mately 13 m/s are observed at 15° N (Figure 4.11). 

By 24 h, the tilt in the sea-level pressure trough has 
increased to produce a stronger pressure gradient in the 
northwest quadrant and a weaker gradient in the northeast 
quadrant of the anticyclone. There has been a corresponding 
change in the wind field (Figure 4.12) . A large area of 
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Figure 4.11. Wind vectors at 12 h for Experiment 
I. 8 mm = 20 m/s. 
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Figure 4.12. Wind vectors at 24 h for Experi- 
ment I. 7 mm = 20 m/s. 
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Figure 4.13. 



Wind vectors at 36 h for Experi- 
ment I. 7 mm = 20 m/s. 
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Figure 4.14. Wind vectors at 48 h for Experi- 
ment I. 8 mm = 20 m/s. 
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Figure 4.15. Wind vectors at 60 h for Experi- 
ment I. 7 mm = 20 m/s. 
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northeasterly winds is present in the southeast quadrant. 
Pressure centers have continued to weaken. A well-defined 
trough has formed along the equator with wind speeds as 
high as 13 m/s. The tilt and amplitude of the potential 
temperature wave have increased. Strong southwesterlies 
ahead of the trough are advecting warm air northward from 
20° N to 45° N. 

At 36 h, the anticyclone has intensified by about 5 mb. 
Ridging dominates the pattern particularly south of 30° N. 

The trough has filled as the flow has become more zonal. 

The cyclone has intensified slightly. Weak cyclonic flow is 
evident at the base of the trough (Figure 4.13). Amplifi- 
cation of the potential temperature field has continued with 
strong warm air advection ahead of the trough. Cold air 
advection behind the trough increased south of 20° N. With 
cold air moving south and warm air moving north, V'0' is 
positive, which indicates a baroclinic process. The area 
and speed of northeasterly winds on the eastern side of the 
anticyclone have increased. 

At 48 h, the ridging south of 30° N has expanded east- 
ward while the anticyclone weakened. North of 45° N the 
cyclone has intensified slightly and enlarged. The pressure 
gradient is oriented north-south. Although not evident on 
the pressure field, a weak cyclone has formed at the base of 
the trough (Figures 4.4 and 4.14). The wave structure in the 
potential temperature has amplified further. The 305° K 
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isotherm has moved to within 3° of the equator. North of 
30° N, the field is 180° out of phase with that south of 
30° N. Winds are northeasterly on the upwind side of the 
thermal trough. 

In the next twelve hours, rapid development has occurred 
in the cyclone. The central pressure has dropped almost 
10 mb. A secondary anticyclone formed 15° northwest of the 
primary high pressure center. The wave in the potential 
temperature field has continued to amplify. The 305° K 
isotherm has reached 1° S. The field does not reflect the 
surface pressure pattern north of 30° N. Northeasterly 
winds of greater than 15 m/s on the eastern side of the 
primary anticyclone are observed in the equatorial region 
(Figure 4.15). The secondary anticyclone is not adequately 
reflected in the wind fields indicating that this feature 
is not in geostrophic balance. 

This experiment demonstrates the effects of a baroclinic 
jet on the stable wave number 4. The jet provided enough 
forcing to produce cold temperatures and strong northerly 
winds in the tropics. The 305° K isotherm moved southward 
from near 10° N to south of the equator. At 60 h, the 
secondary pressure centers indicate growth of wave number 8 
which is unstable. This wave is induced by the nonlinear 
interaction of the long wave with itself. 
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B. EXPERIMENT II 



In this experiment, the model was initialized with the 
baroclinically unstable mean flow and the wave niimber 8 
disturbance described in Chapter III. The model was 
again integrated forward in time to 60 h. Sea-level 
pressure fields for 12, 24, 36 and 48 h are shown in Figures 
4.16 through 4.19. The potential temperature fields for 
12, 24, 36 and 48 h are depicted in Figures 4.20 through 4.23. 

The results of this experiment are used to determine 
when to add the stable wave number 4 flow which will advect 
the synoptic scale disturbance southward in Experiment III. 

One requirement for this simulation is to have a sea-level 
pressure difference of at least 25 mb between high and low 
centers to ensure an adequate northwesterly current for 
surge initiation (Baker, 1983). Another requirement is the 
appearance of a high pressure center representing the 
Siberian anticyclone to the northwest of the low pressure 
center . 

As shown in Figure 4.16, the disturbance is primarily a 
feature of the northern hemisphere mid- and upper-latitudes 
with a very flat pressure gradient in the tropical region. 

In the first twelve hours a low has formed in the trough 
approximately at 60° N. The trough has a slight NE-SW 
tilt and has moved eastward approximately 10° at 45° N. 

A sea-level pressure difference of greater than 25 mb exists 
between high and low centers. The potential temperature and 
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Figure 4.16. Sea-level pressure field at 12 h for 

Experiment II. Contour interval is 5 mb. 
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Figure 4.17. Sea-level pressure field at 24 h for 

Experiment II. Contour interval is 5 mb. 
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Figure 4.18. Sea-level pressure field at 36 h for 

Experiment II. Contour interval is 5 mb. 
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Figure 4.19. Sea-level pressure field at 48 h for 

Experiment II. Contour interval is 5 mb. 
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Figure 4.20. Potential temperature field at 12 h 
Experiment II. Contour interval is 
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Figure 4.21. Potential temperature field at 24 h for 

Experiment II. Contour interval is 5° K. 
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Figure 4.22. Potential temperature field at 36 h for 

Experiment II. Contour interval is 5° K. 
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Figure 4.23. Potential temperature field at 48 h for 

Experiment II. Contour interval is 5° K. 
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wind fields indicate a baroclinic process with warm air 
being advected northward and cold air being advected 
southward. 

At 24 h, the low centers have moved 5° southeast of 
their previous positions and have intensified. A sharp 
surface-pressure ridge extending north to 80° N separates 
the disturbances. Eastward movement of the trough has 
continued at the previous rate. Development of the primary 
disturbance in the wind and potential temperature fields 
parallels that of the pressure field. Some short wave noise 
is apparent in the wind field. Noise is to be expected 
because no surface frictional dissipation, diffusion or 
smoothing techniques are included in the model. These 
techniques tend to dissipate surge energy, and the objective 
of this study is the simulation of. intense monsoon surges. 

From 24 h to 36 h, rapid deepening has occurred in the 
cyclone. The central pressure of the low has dropped at 
least 15 mb. The ridge between low centers has weakened, 
but the sea-level pressure difference between high and low 
centers has increased to more than 50 mb. A broad band of 
low pressure appears to dominate the zone between 50° N and 
60° N. Surface warm and cold fronts are clearly evident 
in the temperature field at 36 h (Figure 4.22). The potential 
temperature shows a poleward flux of warm air and a migration 
of cold air south of 30° N. 
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The rapid deepening of the surface pressure low con- 
tinued between 36 h and 48 h. The central pressure dropped 
approximately 10 mb. High pressure centers formed near 
40° N, and the sea-level pressure difference between the 
high and low pressure centers is greater than 70 mb. The 
northern portion of the strong baroclinic zone is advected 
southeastward about 7° latitude. The potential temperature 
field shows some additional migration of cold air south of 
30° N (Figure 4.23) compared to 12 h earlier. At 60 h, the 
300° K isotherm follows the 30° N latitude circle, and the 
305° K isotherm is a little south of 20° N. 

Analysis of these fields reveals several key features. 

A sea-level pressure difference of greater than 25 mb 
between high and low centers is present from 12 h onward. 

The sea-level pressure difference actually increases 
throughout the run. A high pressure center northwest of the 
low is best depicted at 24 h; it weakens afterwards. A 
strong Siberian high is not observed probably due to the 
presence of a strong jet and lack of topographic and radia- 
tive forcing. The period of most rapid cyclonic develop- 
ment is between 24 h and 36 h. The cold air and strong 
northerly winds never penetrate south of 30° N. 

Based on these results, wave number 4 should be added to 
the growing cyclone between 24 h and 36 h. The most rapid 
cyclogenesis is occurring at this point. Adding the larger 
scale wave will further advect the cold air south to help 
initiate the surge. 



54 



C. EXPERIMENT III 

The initialization for this experiment is the same analy 
tic scheme as Experiment II. However, after 36 h of model 
time, the wave number 4 fields are added. (Note: The 

zonal mean has been removed from the wave number 4 fields.) 
Sea-level pressure, potential temperature, and wind fields 
are depicted in Figures 4.24 through 4.35. 

The development of the disturbance through 36 h is 
described in Experiment II. The effect of wave number 4 on 
the disturbance is readily apparent when comparing the 36 h 
fields of Experiments II and III. A high pressure center 
of at least 1030 mb forms in the ridge. The primary cyclone 
deepens another 35 mb with the addition of wave number 4. 
This rapid deepening produces a pressure. surge . 

The second stage of the monsoon surge is characterized 
by a sharp drop in surface dewpoint (Baker, 1983). This 
model does not include moisture or its effects, but observa- 
tional studies (Chang et , 1979 and Chu, 1978) show that 
this stage can also be associated with a drop in surface 
temperature. At 36 h, a trough of cold air extends south- 
ward to approximately 17° N at the base of the high pressure 
however, the potential temperature gradients are very weak 
south of 30° N. The winds are primarily northerly at 
15 m/s on the eastern side of the anticyclone. 

From 36 h to 48 h the weaker low moved east toward the 
stronger low. The long wave ridge weakened and tilted more 
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Figure 4.24. Sea-level pressure field at 36 h for 
Experiment III. Contour interval is 
10 mb . 
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Figure 4.25. Sea-level pressure field at 48 h for 
Experiment III. Contour interval is 
10 mb . 
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Figure 4.26. Sea-level pressure field at 54 h for 
Experiment III. Contour interval is 
10 mb . 
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Figure 4.27. Sea-level pressure field at 66 h for 
Experiment III. Contour interval is 
10 mb . 
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Figure 4 \ 28 . Potential temperature field at 36 h for 
Experiment III. Contour interval is 
5° K. 
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Figure 4.29. Potential temperature field at 48 h for 
Experiment III. Contour interval is 
5° K. 
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Figure 4.30. Potential temperature field at 54 h for 
Experiment III. Contour interval is 
5" K. 
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Figure 4.31. Potential temperature field at 66 h 
Experiment III. Contour interval is 
5° K. 



'' o 

A' 30 N 




/ 30° S 



/ 



for 



63 



X ' ' ■ '' '' ‘ 7 ' 'v' ' \ 

X/ ^ » - \ V \ V / ^ ^ \ V / V \ 




^ ^ ^ ^ \ V- ^30* S 



. .' f r 

w V, ^ ^ 

\ ^ r f 

' k\V ' / 

6o’s ,v ^",<7 VJm7^ VV_,. .I„3„ 



^ro<wX' 



Figure 4.32. Wind vectors at 36 h for Experiment 
III. 4 nun = 20 m/s. 
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Figure 4.33. Wind vectors at 48 h for Experiment 
III. 8mm=50 m/s. 
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Figure 4.34. Wind vectors at 54 h for Experiment 
III. 8 mm = 50 m/s. 



66 



60 




VV60 N 



90 KJcJga N 

-^>v 

^ v->yy.^/. y yV V 

// N/^i ^ N/ i i 

/*\- / -^>»\^ / \ ^ ^ 




\ k ^ 1 ^ N \ ^ ^ M ^ ^ 

'n , k i ^ 1 I V 



/ ^' 30 S 



V/ 



ik ^ ^ 

go” s V—/ ^ \ > 

a>u> 



V 




'60 S 



Figure 4.35. Wind vectors at 66 h for Experiment 
III. 8 mm = 50 m/s. 
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NE-SW. The small cyclone deepened at least 10 mb. The 
large cyclone deepened approximately 10 mb and moved 
slightly eastward with the mean flow. The tilting and 
weakening of the ridge caused three anticyclonic circulations 
to develop within the ridge particularly toward the north- 
west. At 48 h, winds are northeasterly from 40° N to about 
10° N with a speed of at least 18 m/s. 

The 48 h potential temperature field shows a broad 
tongue of cold air extending southwest near 20° N as if 
spreading over the South China Sea. The furthest extent is 
approximately 15° N. The thermal gradient remains weak south 
of 30° N and especially south of 20° N. A strong baroclinic 
zone is depicted behind the cold front. 

Six hours later the secondary low pressure center is 
located 5° latitude further southeastward. The primary low 
center remained stationary. The ridge elongated and moved 
westward. This movement resulted in a band of northeasterly 
winds (of approximately 20 m/s) extending southwestward 
from 40° N to at least 10° N (Figure 4.34) as if along the 
eastern coast of China. These northeasterly winds advect 
colder temperatures as far south as 15° N (Figure 4.30). 

The strong baroclinic zone moved approximately 2.5° 
southeastward . 

From 54 h to 66 h the low pressure center moves east- 
ward, and the high pressure center moves westward. The 
pressure gradient south of 30° N decreases. The NE-SW 
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band of northeasterly winds narrowed but continues to 
advect cold temperatures as far south as approximately 
10° N. The strong baroclinic zone continued moving south- 
eastward but at a slower rate. 

Several key observations are made based on this analy- 
sis. Northeasterly winds move rapidly south to the tropical 
regions in 12 h. These winds brought cold mid-latitude 
air to the tropics in 30 h. Compared to Experiment I, the 
northeasterly winds are stronger but do not go quite as 
far south into the equatorial region. The sea- level pressure 
pattern progresses eastward very little because additional 
forcing is necessary to overcome the beta effect on wave 
number 4. The pressure ridge and the cold air tongue are 
better defined than the wind surge. Baroclinic processes 
are present in the results of this experiment south of 
30° N. 
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V. CONCLUSIONS 



The interaction of a stable planetary wave with a 
growing cyclonic disturbance in a mid-latitude baroclinic 
flow was examined. The purpose of this study was to deter- 
mine the contribution of a quasi-stationary stable wave 
number 4 to the initiation of a monsoon surge as a wave 
number 8 cyclone develops. 

Several circulation features have been identified as 
characteristic of monsoon surges. They include: (1) the 

strengthening of north-northeasterly winds in the South 
China Sea, (2) a sharp drop in surface temperature at obser- 
vation stations in China and Southeast Asia, (3) the shallow 
structure of the surge, (4) the strong baroclinic zone and 
(5) northerly winds in the tropics. In varying degrees, 
each of these characteristics were observed in one or more 
of the experiments. 

In Experiment I, the jet provided enough additional 
forcing to the wave number 4 to produce cold temperatures 
and strong northerly winds in the tropics. The 305° K 
isotherm moved from near 10° N to south of the equator in 
48 h. The winds were initially northerly with a small area 
of northeasterly winds. Throughout the experiment, the winds 
on the east side of the anticyclone became more northeasterly 
and increased in speed in the equatorial region. A strong 
baroclinic zone was observed on the downwind side of the trough. 
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Experiment II showed the evolution of an unstable baro- 
clinic wave which was limited to the mid-latitudes. The 
300° K isotherm moved south approximately 5° latitude in 
48 h, as opposed to 10° and south of the equator in Experi- 
ment I. A strong baroclinic zone, which is characteristic 
of a frontal zone, was evident in the low-level potential 
temperature fields at the mid-latitudes, but it did not 
propagate southward into the tropics. In this experiment 
there were no significant northerly winds south of 45° N 
and no cold temperatures south of 30° N. 

Several characteristics of monsoon surges were observed 
in Experiment III in which the wave number 4 and wave number 
8 were combined. North-northeasterly winds reached the 
tropical regions 12 h after the addition of the long wave. 
Cold air moved southward to approximately 10° N within 30 h 
after the addition of the long wave. A strong baroclinic 
zone was evident in the potential temperature fields. It 
was advected southward to the tropics. Northerly winds were 
observed in the tropics at least as far south as 10° N. 
Experiment III roughly approximated events occurring during 
a cold surge. 

These results suggest that the combination of quasi- 
stationary stable long wave and a developing mid-latitude 
cyclone initiate a monsoon surge. Additional forcing appears 
to be required to keep wave number 4 stationary. Wave number 
4 also should be modified so that it does not extend to the 
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equator initially- Moving the jet and the maximum lati- 
tudinal amplitude of the perturbation to 35“ N would increase 
the baroclinicity south of 30“ N as well as more closely 
model the observed position of the jet during a surge event. 
This study found that a quasi-stationary stable wave niamber 
4 and an unstable wave number 8 could initiate a monsoon 
surge . 
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APPENDIX A 



MODEL PRIMITIVE EQUAT IONS 



Equation of State: 



pa = RT 



Hydrostatic Equation: 
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Horizontal Momentum Equation: 
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Continuity Equation: 
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1st Law of Thermodynamics: 
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APPENDIX B 



LIST OF SYMBOLS 



A constant 

a radius of the earth 

C constant 

C specific heat at constant pressure 

ir 

f Coriolis parameter 

fQ value of Coriolis parameter at 45° N 

g gravitational acceleration 

m wave number 

p pressure (p = p + p' + P') 

Pq surface pressure at 45° N, 1013.25 mb 
p^ maximum model pressure height, 125 mb 

Pg surface pressure 

p^ pressure at tropopause 

R gas constant for air 

S smoothing operator 

T temperature 

Tq surface temperature at 45° N, 288° K 

Tg surface temperature 

u zonal velocity component (u = u + u' + U') 

Uj^ maximum velocity, surface: 5 m/s, upper: 6 5 m/s 

u^ surface velocity component 
u^ upper level velocity component 
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V 



meridional velocity component (v = v + v' + V) 
Z terrain height 

a specific volume 

3 frequency 

Y halfwidth of jet 

r dry adiabatic lapse rate, 6.5° K/km 
0 potential temperature 

X longitude 

7T 3.14159 

n normalized pressure 

a vertical coordinate 

d vertical velocity in sigma coordinate system 

(j) latitude 

(})q latitude of jet 

$ geopotential, also current terrain height 

angular velocity of the earth 

oj angular velocity of air relative to the earth 
(mean flow) 

ip stream function 
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